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Abstract. We investigate the quantum entanglement dynamics of undriven 
anharmonic (nonlinear) oscillators with quartic potentials. We first consider the 
indirect interaction between two such nonlinear oscillators mediated by a third, linear 
oscillator and show that it leads to a time-varying entanglement of the oscillators, the 
entanglement being strongly influenced by the nonlinear oscillator dynamics. In the 
presence of dissipation, the role of nonlinearity is strongly manifested in the steady 
state dynamics of the indirectly coupled anharmonic oscillators. We further illustrate 
the effect of nonlinearities by studying the coupling between an electromagnetic field 
in a cavity with one movable mirror which is modeled as a nonlinear oscillator. For this 
case we present a full analytical treatment, which is valid in a regime where both the 
nonlinearity and the coupling due to radiation pressure is weak. We show that, without 
the need of any conditional measurements on the cavity field, the state of the movable 
mirror is non-classical as a result of the combined effect of the intrinsic nonlinearity and 
the radiation pressure coupling. This interaction is also shown to be responsible for 
squeezing the movable mirror's position quadrature beyond the minimum uncertainty 
state even when the mirror is initially prepared in its ground state. 
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1. Introduction 

The last decade has seen a surge of interest in investigating the quantum properties 
of micro- or nanomechanical systems [1]. There have been numerous efforts to prepare 
entangled or other non-classical states of such mechanical systems. With the advent of 
techniques such as laser cooling [2] , these mechanical systems can be cooled sufficiently 
close to their ground states and hence tailored to mimic the physics of quantum harmonic 
oscillators to a very good approximation. There are many proposals which aim to 
investigate the dynamics of such mechanical systems operating deep in the quantum 
regime. These include coupling them to other quantum systems such as an ultracold 
atomic Bose-Einstein Condensate (BEC) [21 H], a Cooper pair box [5l |6] and even 
entangling two distant oscillators [5l El [TJ |8] . Often the main motivation is to test the 
foundations of quantum mechanics, but these nano-mechanical systems also turn out 
to be extremely good candidates for applications such as ultra-sensitive measurement 
devices [9]. 

The physics of anharmonic oscillators has been studied in great detail by several 
authors both in the classical and quantum domain [lOl [HI [12l [13] . In particular, Milburn 
has invesigated the quantum and classical dynamics of an anharmonic oscillator in phase 
space [11] and has shown that decoherence induced state reduction results in quantum 
to classical crossover in a nonlinear oscillator. Also in [13] a quantum master equation 
has been derived for a doubly clamped driven nonlinear beam. 

Many of the existing schemes, which aim at exploring the quantum dynamics of 
nano-mechanical systems, treat them as harmonic oscillators. However, an anharmonic 
(nonlinear) oscillator in the quantum regime offers a number of intriguing new 
possibilities for quantum state preparation and manipulation. One of many motivations 
for studying nonlinear oscillators is that by active cooling techniques, such as laser 
cooling, the thermal fluctuations of these nanomechanical systems can only be reduced 
to the standard quantum limit. If a reduction in noise is sought beyond this limit, 
then squeezing the quadratures of these mechanical oscillators is required. For this 
one typically relies on nonlinearities. There already exists many feasible schemes that 
explore the possibility of squeezing the state of a mechanical oscillator [HI [151 [16] . 
Moreover, coherent nonlinear effects are of great interest as they turn out to be important 
resources for processing universal quantum information with continuous variables [T7] . 

In the present work, however, we shall concentrate on the influence of an intrinsic 
nonlinearity on the entangled states of two such indirectly coupled oscillators. We 
address a situation where the anharmonic oscillator is coupled to a second quantum 
system. Firstly we investigate the quantum dynamics of two such anharmonic oscillators 
interacting with a linear oscillator. We show that as a result of indirect interactions 
mediated by the linear oscillator, the two nonlinear oscillators exhibit a time-varying 
entanglement. Interestingly, we find that the effect of nonlinearity is much more 
pronounced for certain initial states. When dissipation is included, the effect of 
nonlinearity strongly governs the steady state evolution of the indirectly coupled 



Quantum entanglement of anharmonic oscillators 



3 



nonlinear oscillators. 

As a second illustration of the effect of the nonlinearities we investigate the unitary 
evolution of a cavity mode interacting with a movable mirror which is modeled as an 
anharmonic oscillator. We provide a full analytical treatment of a physical model that 
describes this interaction in a regime where both the nonlinearity and the coupling due 
to radiation pressure is weak. We show that unitary evolution results in time-dependent 
entanglement between the oscillator and the cavity mode. Moreover, under the joint 
action of radiation pressure coupling and intrinsic nonlinearity the movable mirror will 
also exhibit non-classical dynamics [T8] . 

Nonlinear effects are typically small in nano- cantilevers, since the amplitude of their 
oscillations are inevitably small compared to their length. Moreover, it is difficult to 
control the nonlinearities externally. In this paper we propose to use an electromagnetic 
setup based on a Helmholz coil configuration, where the nonlinearity stems from the 
fact that the energy due to the interaction between the magnetic field produced by the 
coils and permanent magnets at the tips of the cantilevers has a term that depends on 
the fourth power of the defiection of a tip from its equilibrium position. This allows 
us to externally tune the strength of the nonlinearity, which may be difficult in other 
realisations [I5l [16] . 

The paper is organized as follows. In section [2] we introduce the theoretical model 
describing, in detail, the indirect interaction of two anharmonic oscillators mediated 
via a linear oscillator. It is followed by section |3l in which we investigate the coherent 
interaction between an anharmonic oscillator and a quantized cavity mode. In section 
Hlwe briefiy sketch a scheme to induce nonlinearities for the mechanical oscillator, and 
finally we present our conclusions in section |5l 

2. Indirectly coupled anharmonic oscillators 

In this section we shall explore the quantum dynamics of two anharmonic oscillators, 
which interact with the same linear oscillator. We will keep the theoretical treatment 
general at this point, but in section H] we will discuss a potential realisation of the 
required nonlinearites. 

2.1. Unitary Dynamics 

Consider two identical micro- or nanomechanical oscillators each of mass m and 
operating in the quantum regime with fundamental vibrational frequency Um- Denoting 
the position and momentum operators of each oscillator by qi and pi, where i = 1,2, 
the free evolution of the oscillators is governed by the Hamiltonian 



i=l L -I 

If we can modulate the potential seen by the oscillator such that there is an additional 
term proportional to qf, then this will introduce an effective nonlinearity for the 
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mechanical oscillator. The Hamiltonian of two such independent anharmonic oscillators 
then takes the form 



2m 2 



(2) 



where jSqf is the nonlinear interaction energy. Expressing the position and momentum 
operators of each oscillator as 



where a^{a) and 6"'' (6) are the creation (annihilation) operators for the vibron excitations 
of the two anharmonic oscillators, and neglecting all the counter-rotating terms, takes 
the form 

H/h ^ ujmia^a + b^b) + I3{hl + ha) + I3{hl + h^), (3) 

where ha and hi, are the number operators of the two anharmonic oscillators and /3 is 
the nonlinearity strength. It is worth stressing that in general, for a driven nonlinear 
oscillator, the oscillation frequency depends on the driving amplitude [19], although 
a single resonance frequency can still be a valid approximation in the case of a very 
weak driving force. Moreover, in the present work we are considering a system of two 
undriven nonlinear oscillators, for which assigning a single resonance frequency seems to 
be a reasonable assumption. Nonetheless, depending on the initial excitation amplitude, 
the nonlinear oscillator might exhibit multistable behavior. But as long as the initial 
average number of excitations {h) of each oscillator is such that Um + {n)P ~ Um, 
the assumption of a single resonant frequency for each oscillator is still a reasonable 
approximation. Keeping this is mind in the discussion to follow, we shall restrict 
ourselves to low-excitation subspaces of each oscillator. 

We are interested in the indirect interaction between the two nonlinear oscillators 
mediated by a linear oscillator with quantized energy levels equispaced by hjjj. The 
indirect coupling is advantageous because it allows for accurate control of the interaction 
strength by manipulating the mediating oscillator, and consequently gives a handle 
on the quantum dynamics of the two nonlinear oscillators. The importance of the 
indirect interactions can be further appreciated in the dissipative regime. There, if the 
dissipation rate of the mediating oscillator is much faster than the thermal relaxation 
rates of the individual oscillators, then steady state entangled states of the oscillators can 
be achieved. The linear oscillator is here assumed to be addressable by electromagnetic 
radiation created by excess charge or by nano-magnets at the tip of the oscillators, which 
produces an oscillating electromagnetic field [3]. Making the rotating wave and dipole 
approximations, the unitary evolution then corresponds to the Hamiltonian 

Hi/h = ujm{a)a + b^^b) + uc^c + I3{hl^ + ha + hi + hb) 

+ K{a'^c + b^c) + h.c, (4) 
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where c^,c are the creation and annihilation operators for the single quantized mode 
of the linear oscillator, which couples symmetrically — with coupling strength k — 
to each of the nonlinear oscillators. The Hamiltonian (jl]) may, for instance, describe 
the coherent interaction of two anharmonic oscillators with an ultracold Bose-Einstein 
condensate (BEC) [8]. In which case, in the limit of low atomic excitations, the creation 
and annihilation operators and c will be analogous to the collective atomic raising 
and lowering operators J~^, J~ |22l[8]. As shown in |8], the indirect coupling strength k 
between the two nonlinear oscillators can be made to exceed the direct coupling ^direct 
between them. Moreover, the nonlinearity strength /3 can also be made stronger than 
the direct coupling strength ^direct such that ^direct < f3 < k. For instance by following 
the treatment in [8] and treating the two nanocantilevers as anharmonic oscillators with 
a zero-point oscillation amplitude of 50 pm and with a ferromagnet with 10^ atoms 
on each cantilever tip. If the BEC cloud is treated as a linear oscillator trapped at a 
distance d = 1 fim above the nanocantilevers and contains = 10^ atoms in the trap 
center then K/^direct = 8 and (3/ndhect = 5 (see Section 4). 

A general solution of (jl]) may be found, but since the Hamiltonian conserves the 
total number of excitations a significant simplification occurs. If these oscillators can 
be cooled near to their ground states we can restrict ourselves to subspaces with few 
excitations. Furthermore, in order to simplify the analytical and numerical treatment we 
will truncate the Hilbert space of the middle linear oscillator to its two lowest excitation 
subspaces. Previously we have found that this assumption results in a rescaling of the 
Rabi oscillations without altering the qualitative picture [H]. 

In the one-excitation subspace, the relevant basis states for the unitary dynamics 
governed by the Hamiltonian (jl]) are |l)a|0)b|0)c, |0)a|l)b|0)c and |0)a|0)b|l)c. Here 
I l)a|0);,|0)c denotes a state where one of the anharmonic oscillators is in its first 
excited state while the other nonlinear oscillator and the linear oscillator are in their 
ground states, and analogously for the other combinations. If we assume that the 
energy splitting of the linear oscillator can be brought in resonance with the oscillation 
frequency of the anharmonic oscillator, then in the interaction picture the Hamiltonian 
(jl]) takes the form 

Hint = ^P{nl + ha + nl + h^) 

+ hK{a}c + Pc) + h.c. (5) 

General initial states of the nonlinear-linear coupled oscillator system in the subspace 
of one excitation can be written as 

1 

\^lit))= 5^C,,i_,-o(t)|j)a|l-jH|0). 

i=o 

+ Co,0,lW|0)a|0)6|l)c (6) 

In the one-excitation subspace, with the initial condition |\&(0)) = |l)a|0);,|0)c, the time- 
evolved wave function due to the coherent interactions between the two anharmonic 
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oscillators and the effective two-level system becomes 

= «l(t)|l)a|0)b|0), + a2it)\0)a\l)b\0)c + «3W|0)a|0)6|l)e, (7) 



where, 



= (2 + cos(i^it/2) - sm{K,t/2)'j (8) 



«2(t) = ( - 2 + cos(/^it/2) - sin(i^it/2) j (9) 

as{t) = - 2z^e-'^'/'sm{K,t/2)\0)ami)c, (10) 
with i^'i = + 8k^. In the limit /3 — )■ we obtain 

--^=(sin(y2K«))|0)„|0)i|l)„ (11) 

which coincides with the wavefunction that describes the dynamics of two linear 
oscillators interacting with an effective two-level system as discussed in [8]. It should be 
noted that the effect of the nonlinearity cannot be fully appreciated in a one-excitation 
subspace. In this case the effect of the nonlinearity can be mimicked by making the 
two oscillators non-resonant with the mediating two level system. Hence to better 
understand the effect of the intrinsic nonlinearity on the quantum dynamics of each 
oscillator, we have to study also the dynamics of the system in the two- and three- 
excitations subspace. 

With the result of the unitary evolution for all three excitation subspaces in hand, 
we can now attempt to characterize the entanglement between the nano-cantilevers, and 
by doing so try to understand the influence of the inherent nonlinearities. The time- 
dependent state of the two anharmonic oscillators is a mixed state found by tracing 
over the degrees of freedom of the linear oscillator. To quantify the entanglement in a 
bipartite system in an overall mixed state, we use the Peres criterion [23]. We compute 
the negativity for the reduced density matrix, defined as 



n 



A/'=^}^|A.|-ll, (12) 

where |Aj| is the sum of the absolute values of all the eigenvalues of the partially 
transposed reduced density matrix [23] of size n. A non zero value of negativity N" 
ensures inseparability of a bipartite system. 

As a result of the coherent exchange of excitation(s) between the two anharmonic 
oscillators mediated indirectly via the linear oscillator, the two nonlinear oscillators 
become entangled. As shown in figure [1] the system of two nonlinear oscillators exhibit 
time- varying entanglement, and at certain instants the entanglement is found to be 
maximal or nearly maximal in both excitation subspaces. For the sake of comparison 
we also plot the negativity for two indirectly coupled linear oscillators. As can be seen 
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Negativity Negativity 




(a) (b) 

Figure 1. (Color online) Degree of entanglement, as measured by the negativity for 
/3/k = (solid) and /3/k = 0.5 (dashed). The initial states are (a) Ci_o,o(0) = 1 (b) 
Co,o,i(0) = 1. 




Figure 2. (Color online) Degree of entanglement, as measured by the negativity for 
13/k = (solid) and /S/k = 0.5 (dashed). The initial states are (a) C2,o,o(0) = 1 (b) 
Ci,i,o(0) = 1. 




Figure 3. (Color online) Degree of entanglement, as measured by the negativity for 
/3/k = (solid) and /3/k = 0.5 (dashed). The initial states are (a) 6*3,0,0(0) = 1 (b) 
Ci,i,i(0) = 1. 
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Negativity Negativity 




(a) 

Figure 4. (Color online) Degree of entanglement as measured by the negativity for 
/3/k = (solid) and = 0.5 (dashed). The initial states are : in (a) a thermal 
mixture of initial asymmetric states of the three lowest lying excitation subspaces, in 
(b) a thermal mixture of initial symmetric states of the three lowest lying excitation 
subspaces with average thermal occupancy 0.1. 

from figure [1], for an initial state given by |l)a|0)b|0)c, the nonlinear oscillators exhibit 
stronger entanglement as compared to their linear counterparts. On the other hand we 
find that for a symmetric initial state |0)a|0)5|l)c, a stronger nonlinearity strength f3 
leads to a more weakly entangled state of the two oscillators. 

The time evolution of the entanglement for the indirectly coupled nonlinear 
oscillators in subspaces of higher excitations are shown in figure [2] and figure [31 As 
can be seen in subspaces of two and three excitations the effect of a nonlinearity is 
clearly imprinted on the entangled state of the two oscillators. All these results indicate 
that the effect of nonlinearity is much more pronounced for certain initial states. The 
dynamics becoming more complex in subspaces with higher excitations. As mentioned 
before, the particular form of nonlinearity that we are interested in is clearly manifested 
in subspaces with higher excitations. Adding more excitations to the oscillators will 
correspond to approaching the semi classical limit. A more realistic scenario is when the 
oscillators start in a mixed state. As an illustration of this case, we plot the logarithmic 
negativity for an initial mixed state of two indirectly coupled nonlinear oscillators in 
figure m As expected, the degree of entanglement is reduced considerably as compared 
to the case of initially pure states. Furthermore, crucially depending on the initial 
state, a non zero value of /3 may or may not enhance quantum entanglement between 
the oscillators. 

2.2. Dissipation of the oscillator 

Every physical system is susceptible to dissipation. A more realistic approach will 
therefore take decoherence induced by an environment into account. Dissipation can 
either occur through the thermalization of the two anharmonic oscillators or through 
decay in the effective two-level system. As a first approximation we assume that the 
two nonlinear oscillators and the linear oscillator are coupled to independent zero- 
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temperature heat baths with couphng rates 'ja^ and 7c, respectively, the effect of 
dissipation on their evolution — under the Born-Markov approximation — is well 
described by a Lindblad-type master equation of the form 

d — i 

—p=— [Hint, p] + C,aP + /^bP + /^cP- (13) 

Here p is the density matrix of the system, and 

Cap = ^{2apa^ — a)ap — pa} a) (14) 

CbP = ^{2hpP - Php - pPh) (15) 

CcP = ^{2cpc^ — c^cp — pc^c) (16) 

are the Lindblad operators representing the coupling of the two anharmonic oscillators 
and the linear oscillator to their independent zero-temperature heat baths. 

A typical numerical solution of ( !T3|l in the one-excitation subspace is shown in figure 
El As can be seen from the figure, the effect of the nonlinearities is clearly imprinted on 
the entangled state of the two oscillators even when they undergo dissipation. As for 
the case of unitary evolution, the effect of an inherent nonlinearity of the oscillator is 
much more pronounced for certain initial states. 

The intrinsic nonlinearity of the two oscillators has another dramatic effect in the 
sense that it determines the dissipative dynamics of the coupled oscillators in higher 
excitation subspaces. To see this we solve f|T3|) in the two-excitations subspace. This 
time we only allow the mediating linear oscillator to undergo dissipation on the time 
scale of interest. An equivalent problem has been studied in [8], where under similar 
conditions long-lived entangled states of two linear oscillators were achieved. 

If /3 = then a — 6 is a constant of motion of the Hamiltonian ([5]). Exploiting this 
fact one can obtain steady states of the two linear oscillators which are entangled [8]. 
On the other hand, for nonlinear oscillators this does not hold true. Here, depending 
on the initial state, one may or may not see a steady-state entangled state of the two 
nonlinear oscillators develop. 

To prove this statement a numerical solution of f|T3l) is shown in figure [61 As can 
be seen, for an initial asymmetric state the steady state is an entangled state of the two 
nonlinear oscillators while for an initial symmetric state the steady state is separable. 
These observations can thus also be used as an indirect signature of the state of the 
nonlinear oscillator. 

As is clear from this section the degree of inseparability of the two indirectly coupled 
oscillators has a non trivial dependence on both the nonlinearity parameter /3 and the 
initial state. Depending on the initial state a non zero value of nonlinearity strength 
/3 can enhance or suppress the degree of entanglement of indirectly coupled oscillators. 
This interesting behaviour holds true both when the system undergoes unitary evolution 
and in the dissipative regime. 
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Figure 5. (Color online) Degree of entanglement, as measured by the negativity for 
/3/k = (solid) and /3/k = 0.5 (dashed) and 7a, b/^ = 7c/k — 0.1. The initial states 
are (a) Ci,o,o(0) = 1 and (b) Co,o,i(0) = 1. 



Negativity Negativity 




10 20 10 20 



(a) (b) 



Figure 6. (Color online) (a) Degree of entanglement, as measured by the negativity 
for /3/k = (solid) and /3/k ~ 0.5 (dashed), 7c/k = 2 and ja,b/K = 0. The initial 
states are (a) C2,o,o(0) = 1 and (b) Ci_i_o(0) — 1 . 



3. Interaction with a quantized cavity mode 

In the previous section we saw how an intrinsic nonlinearity can strongly affect the 
entanglement between the oscillators. Here we will discuss a second physical scenario 
where nonlinearities are playing a key role in the quantum dynamics. The first setup 
we have in mind is an anharmonic oscillator coupled to a mode of a Fabry-Perot cavity 
with one fixed and one movable mirror. 

In [251 EE] the problem of a cavity with a movable mirror has been discussed in 
great detail. In these studies the movable mirror was treated as a simple harmonic 
oscillator. This work was purely analytical, and showed that the coherent interaction 
of a movable mirror with the cavity mode generates various non-classical states of the 
cavity mode and the mirror. Here we are interested in probing the quantum features of 
an anharmonic oscillator. We shall model the movable mirror as a nonlinear oscillator 
with a nonlinearity proportional to x"^, where x is the displacement of the mirror from 
its equilibrium position. In what follows we study the coherent interaction between a 
single quantized cavity mode and a nonlinear mirror coupled by the radiation pressure 
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We derive a closed analytical expression for the time-evolved state of the cavity 
field and the movable anharmonic mirror, which is valid in the limit of weak nonlinearity 
and low radiation pressure coupling. 

If we assume that leakage of photons through the cavity can be neglected, then the 
main source of decoherence is the coupling of the mirror to its surroundings, which to 
some extent can also be avoided [28]. In what follows we therefore neglect dissipation 
and only consider unitary evolution of the system of of the coupled cavity- and nonlinear- 
mirror system. 

Consider a single quantized cavity mode with creation and annihilation operators 
and k, and resonance frequency Uk = 27rc/L, where L is the length of the cavity. We 
assume that the movable mirror has been cooled near to its ground state and thus is 
operating in its quantum regime. Under the action of cavity photon induced radiation 
pressure, the movable mirror will oscillate about its equilibrium position. If we assume 
that the mirror moves a distance x along the cavity axis such that the displacement 
is much smaller than the wavelength of the cavity mode in one cavity round-trip time, 
then the scattering of photons to other cavity modes can be safely neglected [27]. The 
length of the cavity then becomes L + x so that the resonance frequency of the cavity is 
of the form u'/. = 2'itc/ (L + x). The Hamiltonian of the cavity can then be rewritten as 

Hcav = hu'^k^k = 2'Kh k^k, (17) 

Ij ~\~ X 

which, in a quantum description of the mirror motion, becomes 

Hcav = hhJkk'^k — hgkk^k{a^ + a), (18) 



where it is assumed that x/L ^ 1, and gk = {uik/ L) fij 2mum is the radiation pressure 
coupling constant between the nonlinear mirror and the cavity field. Thus the unitary 
dynamics of the above physical system is governed by the Hamiltonian 

U^jh = Ukk^k + {ujm + P)a)a + P{a^af - guk'^k^a) + a), (19) 

where the nonlinear mirror has been approximated by a quartic anharmonicity as in 
([3]). The Hamiltonian in f fT9|) can be rewritten using the transformation 

i^trans = e^^ac"^, (20) 
where the unitary operator S is given by 

S = ^A;tA;(at_a). (21) 

Consequently the operators a and k transform as 

a^a+ Pk, (22) 

k — >■ ^exp — {a} — a) . (23) 

+ P 
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Neglecting the counter-rotating terms, the transformed Hamiltonian in fl2UI) 
becomes 



To further simphfy the analysis we assume that both the nonlinearity and the radiation- 
pressure coupling are weak, so that quadratic and higher orders terms in g^/ {um + P) can 
be neglected. This can be justified since a cavity of length L ~ 10"'^ m and a movable 
mirror with oscillation frequency ujm ~ 10^ Hz and zero-point oscillation amplitude 50 
pm gives ghl'^m ~ 0.01. Thus (12^ reduces to 



trans 



,2 



?7trans(t) = cxp A;t A; + i^^^^^^t^mf] (26) 



+ (a;^ + /3)ata - ^^^'^ (fctfc)^ + /^(ata)^. (25) 

It should be noted that = k^k and = '^Q, are constants of motion since 
[-f^trans, i^/c] = [-f^trans, ^^a] = 0. The transformed unitary time-evolution operator 
corresponding to i^trans takes the form 

X exp[— ^(aJm + l3)ta^a — i/3t(a^d)^]. 

The corresponding untransformed time evolution operator is U(t) = e'^Utva.nsi't)^^ ■ See 
the Appendix for technical details regarding the exact form of U(t). 

Under the assumption of weak nonlinearity and low radiation pressure coupling, 
U{t) describes the undamped motion of an anharmonic oscillator interacting with a 
cavity mode. If we assume that both the cavity mode and the oscillator are prepared 
in coherent states with amplitudes a and t] respectively, then the state of the combined 
system evolves as 

m)) = m\m) 

oo „ 

X e-'"'^^-V),|r7(t))„ (27) 
where the state of the mirror is a mixture of Fock states given by 

oo 

\m)a = Y.m)re-^'^'^^'-/'e-^^'"^'-y=\m)a, (28) 

with 



n V 

m=0 * 



fl(t) = r^e-'^^^+^)i + n[l - e-^(-'"+/3)*]. (29) 

Ura + P 
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Figure 7. (Color online) (a) Wigner function of the movable mirror initially prepared 
in its ground state and interacting with a cavity mode with(a) f3/{uJm + (3) ~ 10"'* (b) 
P/{uJm + /3) = 0. Initially jap = 1; gk/{uJm + /?) = IQ-^ and (a;„ + ^ ir/A. 



In the limit /3 — )■ we retrieve the result obtained in [25| 126] where the mirror state 
reduces to a mixture of coherent states. It is worth noting that even in the weakly 
nonlinear regime the effect of the nonlinearity is clearly imprinted on the state of the 
movable mirror which is now in a mixture of Fock states. Also evident from ( l28l) is the 
inseparable state of the nonlinear oscillator and the cavity mode. As can be seen from 
( l27|) . the anharmonic oscillator exhibits periodic entanglement with the cavity mode, 
except at certain instants where the state of the oscillator is completely separable from 
the cavity mode. This happens when + P)t = 2g7r for q E N. The reduced density 
matrix of the state of the mirror is thus given by, 

oo 

/5™(t) = Trk(|^(t))(M/(t)|)= e-l"lVl'"^ (30) 

m,q,n=0 



'nl y/ml y/qi 

It is now a straightforward exercise to compute the Wigner function W{X, A*) [29] of 
the mirror, which we plot in figure [71 As can be seen from this figure the negativity of 
the Wigner function W{X, A*) clearly identifies the non-classical state of the mirror. It 
should be contrasted with the case of a linear oscillator interacting with a cavity mode. 
There the state of the mirror is a mixture of coherent states and always characterized 
by a positive Wigner function. The evolution of the mirror into a non-classical state is 
an effect of the combination of an intrinsic nonlinearity of the mirror and the radiation 
pressure coupling with the cavity mode. This feature should be compared with results 
obtained in [26] where it has been shown that only a conditional measurement on the 
cavity mode can project the linear mirror into a non-classical state. One would expect 
the amplitude and the phase quadratures of the movable mirror to be influenced by the 
intrinsic nonlinearity in the mirror. In order to quantify this we define two Hermitian 
operators Q and P, which correspond to the amplitude and phase quadratures of the 
movable mirror and are given by 

Q{t) = {a\t) + a{t)), (32) 
P{t) =i{a\t) -a{t)). (33) 
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0.8 




(a) 



(b) 



Figure 8. (Color online) Time variation of the variance of the mirror quadratures P 
(solid) and Q ( dashed) with the mirror initially prepared in its vacuum state, where 
gk/icurn + /?) - 0.06 and \a\^ = 5. (a) /3/(w„ + /3) = IQ-^ and (b) /3/(w,„ + /?) = 0. 
As a result of coherent interactions with a cavity mode an anharmonic oscillator 
exhibits time dependent squeezing beyond the minimum uncertainty limit in one of 
its quadratures. 



The coherent interaction of the cavity with the anharmonic oscillator should be reflected 
in the variance of the quadratures P and Q defined by 

AP'it) = {P\t)) - (P(t))^ (34) 
AQ\t) = {Q'it)) - m))\ (35) 

We analytically solve for -P^(t) and Q'^{t) and plot the variance of the quadratures P 
and Q in figure [HI As can be seen there, the coherent interaction between a quantized 
cavity mode and an anharmonic oscillator induces a time-dependent squeezing in one 
of the the mirror quadrature beyond the minimum uncertainty limit. 

It is worth pointing out that the squeezing in the variance of the mirror quadratures 
beyond the minimum uncertainty limit is the result of a combined effect of the intrinsic 
nonlinearity and the radiation pressure coupling. This can be understood from the 
fact that if the mirror is initially prepared in its vacuum state then it is known that a 
nonlinearity of the form (I3]) alone cannot induce squeezing in the mirror quadratures [I5] . 
As a result of joint coherent interaction with the cavity mode and intrinsic nonlinearity 
an initial vacuum state of the nonlinear mirror evolves into a mixture of Fock states 
which exhibit time dependent squeezing beyond the minimum uncertainty limit. 

4. Origin of the nonlinearities 

The harmonic oscillator is often the result of an approximation of a more complicated 
potential landscape. Nonlinear force terms are often naturally present in many 
physical systems, but they are of higher order, hence small. Here we shall outline 
a physical scheme for inducing nonlinearity of a mechanical oscillator in the form of 
a nano-cantilever [31], where the nonlinear quartic term appears as a lowest order 
approximation. In order to obtain a nonlinear oscillator, we propose to use a hybrid 
system which relies on the electromagnetic coupling between nano-magnets located at 
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the tip of the cantilever and external magnetic fields. Consider a setup consisting of 
two identical circular magnetic coils of radii R placed a distance oi R/2 apart, with 
their common axis along the x direction. This Helmholtz coil configuration is known to 
produce a very uniform magnetic field near the centre. A nanocantilever cooled near 
to its vibrational ground state and fabricated with a strong ferromagnet of magnetic 
moment fl attached to its tip is placed at the centre of the Helmholtz coil setup. 

The magnetic field experienced by the ferromagnet at the tip of the nanocantilever, 
placed at the centre of the Helmholtz coils, is given by 



B(x) 



/^O'^turns-^ 

2R 



TT — X 



1 



(36) 



where x is the displacement of the tip of the oscillator from the centre and / is the 
current in the pair of coils. Simplifying ( 136|) for x/R << 1 we get 



B{x) 



144 / X N 4 
125 



5RV^ 

Thus the interaction energy of the ferromagnet is given by 



int 



-fi ■ B{x) 



O.SfiQfil f x\ 



R 



r) 



(37) 



(3J 



Representing the quantized motion of the oscillator in terms of creation and annihilation 
operators and a, the interaction Hamiltonian takes the form 



int 



hi3{d 



(39) 



where (3 = 1.28/io/^B-^mag-^ao'^/^-R^! A^mag is the number of atoms in the ferromagnet, 
and ao is the zero point amplitude of the nanomechanical cantilever. 

Using the physical setup described above, a nonlinearity of strength /3 can be 
induced in the nanomechanical oscillator provided the zero point motion of the cantilever 
can be made large (see [32] for a review of the present state of the art manufacturing of 
nanomechanical oscillators). For a set of parameters where i? ~ 80 nm, J ~ 1 mA, 
A'mag ~ 10^, and ao ~ 50 pm, one obtains a nonlinearity strength of the order 
of /3 ~ 250 Hz, where we have neglected any finite size effects stemming form the 
nanomagnet and coils. 



5. Conclusions 



We have studied the dynamics of anharmonic oscillators with a quartic nonlinearity in 
two different physical settings. We have described, in detail, the quantum evolution 
of two such anharmonic oscillators interacting indirectly via an effective two-level 
system. The two-level system could also, for example, be represented by some chosen 
collective excitations of a Bose-Einstein condensate. We have shown that indirect 
coherent interactions cause the two anharmonic oscillators to exhibit time- dependent 
entanglement. Inherent nonlinearities in the nano-mechanical systems are found to 
strongly influence the entangled state of the two oscillators. Interestingly, the effect 
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of nonlinearity is much more pronounced for certain initial states. The signature of 
nonhnearity is clearly imprinted on the entangled state of the two anharmonic oscillators 
even when these oscillators are subject to decoherence. Nonlinearity also plays a 
crucial role in determining the steady state evolution of the indirectly coupled harmonic 
oscillators. 

The coupling strength between two oscillators can be characterized by the 
connectivity |T2]. Connectivity as defined in [12] is the ratio of the coupling strength 
between the oscillators and the frequency difference between them, and diverges in the 
limit of identical oscillators. A high value of the connectivity corresponds to coherent 
exchange of excitations between the oscillators, which is desirable in order to operate in 
the strong coupling regime where coherent interactions supersede all the losses in the 
system. In the particular physical model studied here, however, we have found that a 
larger value of the coupling strength does not always guarantee a strongly entangled 
state of the two oscillators. The strength of these quantum correlations also depends on 
the nonlinearity parameter and the initial state distribution. In addition, a very large 
coupling strength also makes the rotating wave approximation questionable. 

As a second illustration of the effect of nonlinearity, we have studied the coherent 
interaction between a single quantized cavity mode and a weakly nonlinear oscillator in 
the form of a movable mirror. In this case we have been able to find an analytical solution 
for the unitary evolution of the state of such an anharmonic oscillator interacting with 
a single quantized cavity mode. In particular, we have shown that non-classical states 
of the mirror arise as a result of the combination of the radiation pressure coupling 
and the intrinsic nonlinearity in the mirror. A non-classical state of the mirror can 
be generated both for initial ground and coherent states. Unlike in previously studied 
cases, non-classical states of the mirror can be generated without the need of conditional 
measurements [25l[26]. In addition we have shown how squeezing appears in the variance 
of the quadratures beyond the minimum uncertainty state. It should be stressed that 
for an initial ground state of a single nonlinear mirror no squeezing will be generated. 
Squeezing only occurs due to the interaction between the nonlinear mirror and the cavity 
field. 
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Appendix 



6. The evolution operator 



The unitary operator S in (12T|) which is used to transform the Hamiltonian H2 in (l25l) 
gives the corresponding transformed time evolution operator 





Using the Baker-Campbell-Hausdorf expansion [21] together with making the rotating 
wave approximation, and also neglecting quadratic and higher order terms in gc/C, (l4T!l 
simplifies to 
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